Diagnosing space telescope misalignment and jitter using stellar images 
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ABSTRACT 

Accurate knowledge of the telescope's point spread function (PSF) is essential for 
the weak gravitational lensing measurements that hold great promise for cosmological 
constraints. For space telescopes, the PSF may vary with time due to thermal drifts in 
the telescope structure, and/or due to jitter in the spacecraft pointing (ground-based 
telescopes have additional sources of variation). We describe and simulate a procedure 
for using the images of the stars in each exposure to determine the misalignment and 
jitter parameters, and reconstruct the PSF at any point in that exposure's field of 
view. The simulation uses the design of the SNAP^J telescope. Stellar-image data in a 
typical exposure determines secondary-mirror positions as precisely as 20 nm. The PSF 
ellipticities and size, which are the quantities of interest for weak lensing are determined 
to 4.0 x 10" 4 and 2.2 x 10" 4 accuracies respectively in each exposure, sufficient to meet 
weak-lensing requirements. We show that, for the case of a space telescope, the PSF 
estimation errors scale inversely with the square root of the total number of photons 
collected from all the usable stars in the exposure. 

Subject headings: cosmology - gravitational lensing, large-scale structure of the universe 
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1. Introduction 



The accelerated expansion of the universe is one of the most puzzling astrophysical discov- 
ery of the century. The proposed explanations are dark energy, modified gravity and feedback 
from density fluctuations. To explore the mystery, a few large astronomical surveys are under- 
way (DES, ACT0, SPl|) or in planning stages (SNAP, DESTINES, LSST0, EUCLIE@). The 
sensitivity of these surveys to the expansion of the universe comes from both cosmological dis- 
tances and growth of density perturbations as a function of the cosmic time or redshift. The 
probes utilized are Type-la supernova, weak gravitational lensing (WL), baryon acoustic oscil- 
lations, galaxy cluster counting (selected by optical or using Sunyave-Zeldovic effect), and the 
integrated Sachs- Wolfe effect (ISW). Among all these probes, weak lensing is potentially the most 
rewarding one if systematics are well under control. The do minating systematic s for weak lens- 
ing measurements include galaxy shape mea s urement errors (IHuterer et all 120061 ; iHevmans et al 
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This paper is concerned with reducing the systematic errors in galaxy shape measurements. 
For future w eak lensing surveys , the tolerable RMS multip licative calibration error on WL shear is 
about 10 -3 ( Huterer et all 120061 : lAmara &: Refregieii 120071 ) . Additive errors in galaxy shear should 
also be held to < 10 -3 ' 5 . Mis-estimation of the PSF will propagate into systematic errors in the 
shear. The size of the PSF must therefore be determined to better than 1 part in 10 3 to avoid 
unacceptable multiplicative shear error; likewise the PSF ellipticity must be know n to 10~ 3 
better to avoid unacceptable additive shear systematic (jPaulin-Henriksson et alll2007l ). 



or 



Unfortunately the PSF of real telescopes changes with time, as well as with color and field 
position. Effects that could change the PSF include thermal expansion of mirrors and support 
structures; pointing jitter due to structural vibrations and tracking errors; for ground-based tele- 
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scopes there are additionally gravity loading, atmospheric distortions, and wind loading. Careful 
engineering of the telescope, mount, and housing can minimize these effects, but there are limits to 
what can be achieved with finite resources. Even for a space telescope it is prohibitively expensive 
to guarantee PSF stability to < 1 part in 10 3 . On the other hand there will be stars in each 
image to diagnose the PSF behavior during each exposure. Recall that the WL analysis requires 
part-per-thousand knowledge of the PSF at each location and each exposure, not necessarily that 
the PSF be constant to this level. In this work, we study how well the PSF can be constrained 
using stellar images, using the proposed space-based SNAP telescope study. 

The outline of the paper is as follows. We present an overview of the task we are trying 
to accomplish in £}2j In £}3j we provide details of the modeling of the SNAP PSF. We describe 
the fitting procedure for the misalignment and jitter parameters in §4l The results of the fit are 
presented in Sj5] and we conclude in £|6l 



2. Stellar "Morphometry" 



Weak lensing measurements aim to extract a map of the cosmic shear from the coheren t 



distortions in the shapes of many distant galaxies ( Kaiser 1998 : Bartelmann Sz Schneider] boOll ). 



Observed galaxy shapes are distorted by the telescope PSF. To lowest order, this can be corrected 
if the PSF across the telescope field of view is known. The PSF can be inferred from the observed 
shapes of foreground stars of suitable magnitude. Because the PSF will drift over time, it is desirable 
to measure the PSF across the field of view for every exposure, using the stars that are interspersed 
throughout the image along with the distant galaxies of interest. We refer to this procedure as 
" stellar morphometry" . 



2.1. Quantities of Interest 

The requirements for a weak-lensing survey can be most simply stated as limits on the tolerable 
error in the second moments of the PSF. We will measure the Gaussian-weighted moments defined 
as the zeroth moment 

M = J dxdyPSF(x,y)W(x,y), (1) 

Jf f dxdyx PSF(x, y)W(x, y) , (2) 
y = j^Jdxdyy PSF(x, y)W(x, y) , (3) 



the first moments 



and the second moments 



Pxx = wJ dx dy (x " ^ 2p5i?(x ' y) w{x > y) ' (4) 
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Pyy = w J dx dy {y ~ y )2pSF{x > y)w{x > y) ' (5) 

Pxy = W Q j dxd y ( x - *)(y - y) pSF ( x iy) w ( x iv) ■ ( 6 ) 



These are calculated under a weight function which we take to be 



■morn 



The width of the weighting Gaussian is chosen to be two times the Airy radius, 

Pmom = 2 X I- 22 — , (8) 

where A is the wavelength of the incident light, / is the focal length and D is the telescope aperture. 
In analogy, we can write down the third moments of the PSF, P XX x, Pxxy, Pxyy, and Pyyy 

We also compute quantities derived from the second moments. The ellipticities and stellar size 

Pxx Pyy 2P X y 2 Pxx ~\~ Pyy 

e ^ = P +P > e * = P +P ;°£ = — y^- ( 9 ) 

1 xx \ 1 yy 1 xx y 1 yy * 

These quantities appear in many approaches to weak-lensing shear measurement Heymans et all 
(|2006l ) and references therein]. The true PSF of an exposure will depend on the field position of 
the star (cc*,y*) yielding ei(x*,y*), e2(x*, y*), <t*(x*, y*). Our goal is to produce an accurate model 
estimate ei(a;*, y*), etc. We will evaluate our success by calculating the RMS residual errors in the 
ellipticity models, 

( e RM S) 2 _ ^ [ei(x ^) _ e^y*)] 2 ) , (10) 
( e RM S) 2 ^ ^ [e2(x ^) _ e 2 (x„y,)] 2 ) , (11) 
and the fractional residual error in the PSF size, 

(^ RMS ) 2 - ^ ^ '- (12) 



2.2. Parametric Models 

If the physical state of the telescope can be described by a small number of time-variable pa- 
rameters {pi}, then the PSF is some function PSF(x,y\pi,x in y ir ) of focal-plane position, telescope 
state and the field position of the star. We note that a great advantage of space-based observatories 
for WL work is that the stability of their environment allows us to engineer a telescope for which 
only a small number of degrees of freedom will vary significantly. For the SNAP telescope, the 
engineering specifications are that all optical systems are stable to well below the WL specification, 
except for: 
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• The alignment of the secondary mirror, which may vary due to imperfect performance of the 
feedback system that stabilizes the temperature of the mirror support structure; 

• The telescope line of sight (LOS) may vary during an exposure and smear the PSF due to 
noise and the finite bandwidth of the attitude control system (ACS) or due to high-frequency 
reaction wheel vibrations that transfer through the structure to optical elements, particularly 
the secondary mirror. 



In ^21 we describe in detail the model we adopt for these disturbances. 



Ground-based telescopes pose a more difficult challenge for PSF modeling, because they have 
a very large number of time-varying degrees of freedom. Indeed the atmospheric distortions have 
infinite degrees of f reedom, forma l ly. O ur analysis thus cannot be considered valid for ground- 
based observatories. IJarvis &: Jainl ([2004]) propose instead that a principal-components analysis be 
performed on the ensemble of PSF patterns observed by the telescope, so that the coe fficients of 
some finite number of principal components become the parameters for the PSF model. I Jain et al 
(J2006|) discuss how changes in Seidel aberrations would be manifested as PSF-change patterns, and 
might serve as a parameter set for PSF modeling. The success of these methods will depend upon 
how well-behaved the telescope and atmosphere are, particularly whether the optically significant 
perturbat i ons a re described by a small number of variables. An alternative approach, applied by 
Wittmanl (|2005l ) to PSF ellipticities induced by the atmosphere, is to determine by a priori analysis 
that the disturbance will be below the WL threshold. 



2.3. Simulations 

We simulate the following strategy: 

1. Locate the stellar images in each exposure. 

2. Measure PSF quantities at each star location; in our case, the second and third Gaussian- 
weighted moments. 

3. Find the PSF parameters {pi} that best reproduce the stellar data. 

4. Use the {pi} and the PSF model to derive the desired ex, etc., at any location in the focal 
plane. 

In the simulation we can then evaluate the RMS residual errors of the PSF model. 

In the simple case where the PSF does not depend on the field position of stars, we can 
consider this procedure to be, essentially, averaging the measured PSFs (and moments) of the 
observed stars. In this case, we expected the error on the PSF moments to be determined by the 
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(quadrature) sum of the signal-to-noise (S/N) levels of all the available stars. In particular, if all 
the stars are dominated by Poisson noise from source photons, then S/N = y/N~, where iV 7 is the 
total number of photons collected from all usable stars in the exposure. We can therefore expect 
that 

e RMS „ g RMS ^ ^RMS = ^-1/2 ^ (13) 

where a is some coefficient of order unity (a = l/v2 for a Gaussian PSF). 

More realistically, the PSF does depend on the field position of stars. In this case, we are 
using the PSF parametric model as a means of interpolating between the stellar positions. If the 
parameters are not too numerous, and cause readily distinguished PSF patterns on the focal plane, 
then we expect equationfTBl to continue to hold. Our simulation will show that this is indeed the 
case for the SNAP design, and we will aim to estimate the coefficient a. 



Paulin-Henriksson et all (|2007l ) instead consider the PSF to be locally constant, and ask how 
large a region will contain enough stars to adequately constrain this locally-constant PSF. They 
then consider this region size to be smallest on which WL observations can be successful. In reality 
both the time-varying PSF contamination and the WL signal will have non-trivial angular power 



spectr a and we have to compare the bandpowers of each in setting our specifications. IStabenau et al 



(|2007l ) investigate how PSF time variation will translate into multipole patterns for a SNAP-like 
sky-scan strategy. In this paper we will simply calculate the RMS residual PSF errors, and note 
that they have characteristic angular scales similar to the telescope field of view. 



3. Modeling the PSF 

3.1. Optical Aberrations and Diffraction 

We adopt standard scalar diffraction theory to evaluate the optical contribution to the PSF. 
The wavefront on the focal plane {x, y} generated by a point source is, 

U(x, y )=c f f d^d v P(C, rjyk-oPD(i:,ff) e -ik(e^+ey V ) e ik^+y V yf _ ^ 



Here £ is the coordinate on the entrance pupil with components £ and r/, C is an uninteresting 
constant, P(£,r]) is the entrance pupil function, k = 2tt/\ where A is the wavelength of the band- 
limited optical light used to form the image, and / is the effective focal length of the telescope optics. 
OPD(£; 6) is the opti cal path di fference caused by the lens/mirrors system which we expand using 



Zernike polynomials (|NoHlll976l ) as basis. The second exponential describes the phase differences 
caused by the off-axis incident light ray in the direction 9 and the third exponential is the phase 
differences caused by the different distance light has to travel beyond the lens/mirrors and reach 
the focal plane. The optical point spread function is 

PSF(x) = \U(x,y)\ 2 . (15) 
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Figure[L] shows the pupil function of SNAP telescope. 
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Fig. 1. — Pupil function of SNAP telescope. The outer radius is 1 meter and the inner radius is 
0.35 meter. The width of the three secondary mirror supporting struts is 4 cm. 



3.2. Optical Path Difference (OPD) 

The OPD map of the perfectly aligned telescope and the derivatives with respect to misalign- 
ment parameters are calculated using ray tracing through the telescope's optical system. The OPD 
is projected onto a Zernike basis, with results shown in Table[TJ 



OPD((, 



n=l 



C n 0) + Y,D m C' mn 0) 



(16) 



where Z n is the n th Zernike polynomial, D m is the the m th misalignment parameter, C n {9) is the 
OPD's Zernike coefficients for pristine telescope and C mn (9) is the Zernike coefficients of do> l> 



dD„ 



In the SNAP telescope design (|Lamptonll2002l ; ISholl et alll2004l ). the secondary mirror position 
is expected to be the only optical dimension to vary significantly with time due to thermal drift. 
The secondary mirror has 5 degrees of freedom (DOF) which include shifts D x and D y in transverse 
directions, the defocus D z , plus tilts around the x axis (T x ) and y axis (T y ). 

In practice we find that C' mn {6) is a very weak function of field for these parameters, and hence 
for the current simulation we take C to be constants. The pristine-telescope Zernike coefficients 
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C n {9) retain field dependence, but the axisymmetry of telescope design reduces the freedom to the 
radial direction. 



Table 1: Zernike coefficients of the OPD in nm for a perfectly aligned telescope with field location 
10.4 mrad off axis C n , and the derivatives C' mn (see equation fTo]) . 



Zernike 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


C n 


7.75 


0.00 


10.87 


0.00 


8.43 


0.00 


-4.45 


0.00 


9.15 


-0.88 


Ac (A*m) 


-2.54 


0.00 


0.14 


0.00 


0.18 


0.00 


1.86 


0.00 


0.00 


0.00 


D y (/an) 


0.00 


-2.55 


0.00 


-0.18 


0.00 


1.86 


0.00 


0.00 


0.00 


0.00 


D z (fim) 


-0.08 


0.00 


-23.68 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.24 


T x (//rad) 


0.00 


1.27 


0.00 


0.63 


0.00 


-0.93 


0.00 


0.00 


0.00 


0.00 


T y (/irad) 


-1.27 


0.00 


-0.48 


0.00 


0.64 


0.00 


0.93 


0.00 


0.00 


0.00 



3.3. Charge Diffusion 

The optical PSF must be convolved with the charge-diffusion pattern of the CCD detector. 
Charge diffusion is modeled as a Gaussian with fixed charge diffusion length = 4/im. If the 
charge-diffusion length were free to vary, it would be degenerate with an isotropic telescope jitter 
(see below). The PSF after charge diffusion is 

PSF(x,y) = [ dx'dy>PSF (x>, y')^^e W {- {x - X ' ? ^ V - ^ ) , (17) 

where PSFq is the optical PSF. We execute the convolution with Fast Fourier Transforms (FFTs). 

3.4. Jitter 

Guiding errors and mirror vibrations, also known as jitter, alter the effective PSF of a finite- 
length exposure. Each exposure hence has a unique PSF map, even if the optics are otherwise 
stable. If the observatory is free to rotate on all three for a space-borne observatory or 

an alt-az terrestrial telescope, then the effect of jitter varies across the field of view, and is not a 
simple convolution of the image with a fixed kernel. Stellar images in the exposure can be used to 
infer the full field dependence of the jitter on the PSF. We demonstrate here that as few as two 
stars are sufficient to fully reconstruct the jittered PSFs, as long as the jitter amplitude is much 
less than the width of the PSF. 
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3.4.1. Effect of Jitter on the PSF 



Assume that the modulation transfer function (MTF, the Fourier transform of the PSF) at 
x = (x,y) from the optic axis is known to be To(k) in the absence of telescope jitter. If the jitter 
has displaced the stellar image by some amount Ax = (Ax, Ay) — which varies in time — then the 
transfer function becomes 

T(k,t) = T (k)e^ k - AxW . (18) 
The PSF for stellar images in the integrated exposure is the time-averaged value 

T(k) = ±£ dtT(k, t) = T (k) <V* k - Ax W^> . (19) 

If k • Ax <C 1, then the exponential can be approximated by a Taylor expansion: 

T(k) = T (k) [1 - ik ■ (Ax(t)) - k T <AxAx T > k/2] . (20) 



The effect of the jitter on the PSF is then fully described by the mean displacement (Ax) 
and by the covariance matrix Cax = (AxAx T ). Further detail of the jitter history is irrelevant. 
The linear term is simply a displacement of the entire PSF, and the quadratic term describes a 
convolution of the jitter-free PSF with a very narrow jitter kernel. With a telescope of diameter 
D, focal length /, and wavelength A, physical optics forces MTF=0 for k < 2D/(Xf). Therefore 
the Taylor expansion is valid if 

in other words the jitter must be much less than the size of the Airy disk. 



3.4.2. Field Dependence of the Jitter MTF 

If the observatory axis is displaced by angles 6 = (9 x ,9 y ,9 z ), i.e. pitch, yaw, and roll, then 
the displacement of the image of star i at (xj,yj) is 

Ax, = f9 x - Vi 9 z (22) 

Ay, = f9 y + Xi 9 z . (23) 

If the roll is non-zero, then the image displacement is field-dependent. Now considering the 
observatory misalignment (jitter) to be a function of time, 

(Ax^Ax;) = f 2 (6 X 6 X ) - 2fy t (9 X 9 Z ) + y\ (9 Z 9 Z ) (24) 

(AxiAy,) = f 2 (9 x 9 y ) + f Xi (9 X 9 Z ) - fy t (9 y 9 z ) - x t y t (9 Z 9 Z ) (25) 
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(A yi A yi ) = f (0 y 9 y ) + 2f Xi (6 y 9 z ) + x\ (9 Z 9 Z ) 



(26) 



In the small-jitter limit, therefore, we find that the effect of the jitter on the PSF at every 
point in the field of view is fully described by the six independent elements of the jitter covariance 
matrix C e = (60 T ). We can write 

T(k, x) = T (k, x) [1 - k • C Ax (C e , x) • k/2] . (27) 

In practice, therefore, if we have an exposure for which the jitter-free PSF pattern is well deter- 
mined, then we can completely determine the jittered PSF anywhere in the focal plane by knowing 
the elements of Cg. The jitter covariance matrix could be determined from perfect knowledge of 
the PSF of any two stars in the exposure. If there are a larger number of stars in the exposure, 
then the six jitter covariances are highly over determined, hence they are easily derived for every 
exposure, even if there are other degrees of freedom in To which must be determined from these 
stars. 

If the jitter amplitude is not small, then there can be a much larger number of moments of the 
jitter history that are important, and a finite number of stellar images may not in general recover 
full knowledge of the effect of jitter on the PSF over the field. We will defer consideration of this 
limit for another paper. 

4. Simulation Procedure 

The simulation process is to: assume fiducial PSF parameters (misalignment and/or jitter); 
create simulated stellar images across the instrumented field of view; measure moments of these 
images; fit the PSF model to these moments; and finally, evaluate the quality of the fit. 

4.1. Fiducial Model 

We analyze a fiducial case in which the secondary mirror is translated 1 fim and rotated 1 
firad from its correct position. This would be considered a very large error for the optomechanical 
system. We have verified that the choice of fiducial model does not influence the RMS residuals to 
the fit. 

When analyzing jitter, we assume an RMS motion of 36 mas in pitch and yaw with 700 mas 
RMS in roll (as expected from SNAP telescope design). The RMS cross-correlation between axes 
is taken to be small by comparison. Again, these do not strongly affect the results. 
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4.2. Simulated Stars 



We use the star counts from the COSMOS H ST survey as rep resentative of high galactic 
latitude fields. The star counts in the COSMOS field (jRobin et afeOQq ) are well fit by (see Figure[2|) 



dn* 
dm 



18.5 x l0 a089m deg-Vag" 1 , 



(28) 



where m is stellar magnitude in the F814W band. 



3000 



:osmos + 
At 



s 

s 



1000 
900 
800 



,„.,,. , „0.08930 ra -1 , -5 

18.4715 x 10 mag deg 



' + + 

+ +S 



18 19 20 21 22 23 24 



Fig. 2. — COSMOS survey star magnitude distribution (jRobin et all l2006) 



If the star is too bright, it would saturate the SNAP CCDs. We hence conservatively assume 
that only stars with 19 < m < 23 will be used for morphometry. The bright limit roughly corre- 
sponds to 50,000 photons for a 300 second exposure. The faint limit corresponds to « 1200 photons 
per star, which is comparable to the sky background in the exposure. Fainter stars will contribute 
little to the PSF knowledge. 

The instrumented SNAP focal plane area is ~ 0.7 deg 2 . This places approximately 2100 mea- 
surable stars on the focal plane, with a total photon count of _/V 7 ~ 2.3 x 10 7 per 300 sec exposure. 
This suggests that an ideal morphometry process would yield e RMS , V2a^ MS ~ 2 x 10 4 , well below 
the required weak lensing specification as discussed in §[H 

We generate a sample of 2100 stars with random positions across the focal plane. The mag- 
nitude of the stars is randomly generated according to the magnitude distribution of equationl281 
The PSF, including optical distortions and charge diffusion is computed, and is used to determine 
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the mean number of photons detected in the CCD pixels. The number of detected photons per pixel 
is drawn from a Poisson distribution, and the quantum efficiency and gain of the SNAP pixels are 
used to compute pixel values. Random dark noise of 5 photo-electrons (as expected in the SNAP 
CCDs) is added to each pixel. The PSF moments defined in equationsQ3[6] and equation[9] are then 
computed from the pixel values. 

4.3. Fitting Process 

For each available star in an exposure, we calculate the second and third moments of a PSF 
image to which shot noise, sky noise, and read noise have been added. With the resulting PSF 
moments as data, we search for the best fit misalignment parameters by minimizing x 2 

Nmom , ^rjmom J .'Yjmoms2 

x 2 =e e ( ■ y ■ ] ■ <»> 

i=l jmorn=l moment 

where N* is the total number of stars, Nmom is the number of independent PSF moments per star 
(Nmom = 3 if only 2nd moments are used), M is star moment (with noise), M is star moment 
calculated from the model (no noise), and o~ moment is 

the rms of star moments. In general, o~ m oment 
depends on star magnitude, filter band and the position of the star on the focal plane. We assume 
a fixed filter band and neglect the (< 10%) dependence on position. We produce a lookup table of 
^moment vs source magnitude by Monte Carlo methods before doing the fit. 

The PSF moments depend nonlinearly on misalignment parameters. As an example, Figure[3] 
shows the dependence of P xx on the defocus parameter D z . The model-fitting procedure is hence 
nonlinear, so slower than a linear x 2 fit- 

5. Results 

In this section, we show results of simulating and fitting for secondary mirror misalignment 
only, and simulating and fitting with the jitter parameters jointly. We also show that the inclusion 
of third moments of the PSFs improves the fit. 

5.1. Fit for Secondary Misalignments Only 

5.1.1. Using Second Moments of PSFs 

We first simulate a single exposure of SNAP, using only the PSF second moments to constrain 
the secondary-mirror misalignment. The fiducial misalignment parameters, their fitted values, 
and the 1 — a uncertainties are listed in Tabled D z is precise to < 10 nm which is well below 
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Fig. 3. — Dependence of PSF moment P xx on z-misalignment (in meters). The star is located on 
the diagonal (i.e. x = y). 



the achievable mechanical stability. Hence the morphometric information can greatly improve 
knowledge of defocus. The other parameters are as much as 20 times less precise: accuracies of 
~ 0.2 /jm in D x and D y seem quite disappointing, for example. 

To better understand this wide range of precisions, we examine the eigenvectors and eigenvalues 
of the parameter covariance matrix, shown in TableEJ The best-constrained eigenvector is along 
the D z direction, i.e. defocus. Two eigenvectors are ~ 80 times more poorly constrained, but this 
is not hard to understand. If the secondary mirror were spherical with radius R, it would have only 
three degrees of freedom, as misalignments with D x = —RT y and D y = —RT X leave the sphere 
invariant. A non-spherical secondary mirror breaks these degeneracies and results in finite (but 
poor) constraints on these eigenvectors. 

The "Residual" row in the table gives the projection of the fitting error onto the eigenvectors. 
It shows that the fitted parameters deviate from the fiducial mostly in the directions that have 
weak constraints. All the deviations are at or below 1.6 a level which is a sign that the fitter is 
working properly. 

The quantity that matters to weak lensing is the precision of the PSF knowledge. The appar- 
ent poor fit to some of the misalignments reflects the insensitivity of the PSF moments to certain 
combinations. This insensitivity also means, however, that poor knowledge of these eigenvectors 
has little adverse effect on our PSF model. We compare the noiseless PSF ellipticities and size 
generated using the fitted misalignment parameters with these generated using fiducial misalign- 
ment parameters. Figured] (left panel) shows the distributions of the residual PSF moments for a 
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Table 2: Misalignment fitting results 



misalignment 


fiducial value 


fitted value a 


1 — a error b 


fit incl. 3rd moments 


1 — a error d 


D x (/im) 


1 


0.8542 


0.1839 


1.0052 


0.1494 


D y (/im) 


-1 


-1.3869 


0.2020 


-1.0719 


0.1602 


D z (/im) 


1 


0.9870 


0.0086 


1.0005 


0.0083 


T x (/trad) 


1 


0.4500 


0.3280 


0.8104 


0.3125 


T y (/trad) 


1 


1.2370 


0.3063 


0.8549 


0.2937 



"Fit using second moments of the PSFs. 

b l-a error of the fit using second moments of the PSFs. 

c Fit using both second and third moments of the PSFs. 

d l-a error of the fit using both second and third moments of the PSFs. 



Table 3: Eigen values and eigen vectors of the misalignment covariance matrix 





eigen 1 


eigen 2 


eigen 3 


eigen 4 


eigen 5 


D x (/im) 


0.0661 


0.4852 


0.7423 


0.4573 


0.0037 


D y (/im) 


-0.5022 


0.0535 


0.4597 


-0.7304 


-0.0107 


D z (/im) 


0.0033 


0.0189 


-0.0058 


-0.0191 


0.9996 


T x (/trad) 


-0.8502 


0.1516 


-0.2812 


0.4185 


0.0063 


T y (/trad) 


-0.1438 


-0.8593 


0.3981 


0.2862 


0.0245 


eigen values (Aj x 10 6 ) 


0.3781 


0.3479 


0.0797 


0.0758 


0.0052 


Residual (xlO 6 ) a 


0.6180 


-0.3789 


-0.0370 


0.0538 


-0.0071 



"Projection of the vector that points from the fiducial parameters to the fitted parameters onto the eigen vectors. 

realization of a single exposure. We see that the PSF errors are well below the 10~ 3 target for all 
three quantities. The spread in PSF errors reflects the variation across the focal plane. 

5.1.2. Including Higher Moments of PSFs 

There is information in the higher moments of the PSF. This could further constrain the 
telescope misalignments. As shown in Tabled the fitted misalignments are noticeably closer to the 
true values and the error bars are reduced when the observed PSF second and third moments are 
used in the misalignment fit. Figure^] shows the reduction in PSF ellipticity and size residuals. 
With the inclusion of third moments in the fit, the reduction of the residual moments as shown 
in Figure^] are much more impressive than that of the one sigma errors shown in Tabled This is 
the manifestation that the contributions to the one sigma errors are dominated by the two near 
degenerate eigens which have little or no effect on the residual PSF moments. In the following, we 
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Fig. 4. — Distributions of the residual PSF moments as defined in the text. Left panel: the fit is 
done using 2nd moments of the PSFs. Right panel: both 2nd and 3rd moments of the PSFs are 
used. 

include third moments in the fit unless stated otherwise. 



5.2. Fitting Misalignment and Jitter Parameters Jointly 

Tabled] shows the results of fitting secondary misalignment and jitter parameters jointly. 
Adding the 6 jitter degrees of freedom to the model roughly doubles the uncertainties on the 5 
misalignment degrees of freedom. The parameter D z is still very precise (at 20 nm level). The 
jitter parameters are determined without significant degeneracies among each other nor with the 
misalignments; this is because these parameters influence the PSF moments with rather distinct 
dependences on field angle. 

The residual PSF ellipticities and sizes for one realization of joint misalignment /jitter fitting are 
shown in Figure(5j As mentioned before, the spreads of the residual PSF e%, e2, and <r* distributions 
in this single realization are due to field dependence. Different realizations of the data (star locations 
and random noise change) produce different mean and spread of the residual moments distribution. 
Figure[6] left panel shows the distributions of e R and cr^ from 50 realizations. We produce a 
single measure of the efficacy of the morphometry procedure by averaging the RMS PSF residuals 
over all focal-plane positions in many realizations, as per equations [T0lll2l The mean values (in 
quadrature) of the distributions in Figured] left panel are exactly those. They are labeled by the 
arrows in the plot and listed in TableE] as well. Calculated using equation ll3l the a values are also 
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Table 4: Fitting for misalignment and jitter parameters 



parameters 


fit incl. 3rd moments 


fiducial values 


1 — u error 


D x (M m ) 


U.OOO ( 


i 
i 


U.ooyo 


D y (fim) 


n i oqo 
-U-laoz 


-i 


yj.oiuo 


D z (//m) 


u.y i zo 


i 
i 


u.uzuo 


T x (>rad) 


2.6843 


1 


0.7503 


T y (//rad) 


1.9457 


1 


0.7207 


(6 X X ) (firad 2 ) 


0.0302 


0.0300 


0.0002 


(9 x 9 y ) (jurad 2 ) 


0.0013 


0.0010 


0.0001 


(6 x e z ) (/irad 2 ) 


0.0049 


0.0010 


0.0045 


{9 y e y ) (/irad 2 ) 


0.0298 


0.0300 


0.0002 


(9 y 9 z ) ( M rad 2 ) 


-0.0011 


0.0010 


0.0046 


(9 Z 9 Z ) (^rad 2 ) 


19.5098 


20.0000 


0.9244 



tabulated in Table[5j Taking the average value, we find a ~ 1.8. So we have 

e RMS ^ e RMS ^ ^RMS = Lg ^-l/3 (30) 

The ensemble average residuals are consistent with zero, i.e. the PSF models are unbiased, to the 
2.5 x 10~ 5 accuracy of our 50 realizations. 



Table 5: Average e and cr^ from 50 realizations and a values 





2100 stars per realization 


A^ 7 = 2.37 x 10 7 


100 stars per realization 


iV 7 = 1.29 x 10 6 




average 


a 


average 


a 


C RMS 

C RMS 

_RMS 
u + 


4.0 x 10" 4 
3.9 x 10" 4 
2.2 x 10" 4 


1.95 
1.90 
1.52 


1.8 x 10~ 3 
1.7 x 10" 3 
1.0 x 10~ 3 


2.04 
1.93 
1.14 



To test our hypothesis that the PSF errors will scale as iV 7 , we repeat the simulation using 
100 stars instead of the fiducial 2100 stars. The distributions of e RMS and (tJ" 1 ^ are shown in 
the right panel of FigureO Again, the quadrature means of these distributions are labeled by the 
arrows in the plot and listed in Table[5j We find a ~ 1.7. So the RMS residuals do indeed scale as 
expected. 

From Figure [5] it is clear that part of the residual errors in the PSF ellipticities are from a 
shift in the mean across the field of view, while the rest is from errors that vary across the field of 
view. The two different types of PSF modeling errors will propagate into different angular scales 
in the WL power spectrum. For the SNAP simulation, we find that roughly half of the e modeling 
variance is in the mean across the field of view. 

Since essentially all the residual variance arises from shot noise in the stars, it will be uncor- 
rected from exposure to exposure. 
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Fig. 5. — Distributions of the residual PSF moments for fit including 5 secondary mirror misalign- 
ment and 6 jitter parameters. Both 2nd and 3rd moments of the PSF are utilized. 

6. Conclusion and Discussion 

Our simulation of "morphometry" for the SNAP telescope demonstrates that the ~ 2000 well- 
measured stars in a typical exposure contain sufficient information to reduce the errors in the 
modeled PSF ellipticity and size to 4.0 x 10 and 2.2 x 10 -4 respectively, giving significant margin 
to meet the ~ 10~ 3 level needed to reduce weak-lensing systematic errors below statistical errors 
of future surveys. For the SNAP telescope design and focal plane, we find e R ~ 1.8/ and 

PSF estimation error in morphometry will be only part of the shape-measurement error bud- 
get, so this margin is important. Other potential source of errors in PSF estimation include charge 
transfer inefficiency (CTI), data compression artifacts, and chromatic PSF dependence that causes 
galaxies' PSFs to differ from stellar PSFs. Shape measurement errors can also arise in the deconvo- 



lution process even if the PSF is known precisely (IHeymans et alll2006l : lMassey et alil2007l ). To have 



a successful weak lensing mission, the sum of these errors must meet the weak lensing requirement. 

We have also assumed that the only time-variable aspects of the PSF are the secondary mirror 
alignment, and small pointing jitter. The SNAP spacecraft is engineered to take advantage of 
the extremely stable space environment so that these are the only relevant degrees of freedom. A 
ground based observatory would suffer in addition the effects of wind, gravity loading, and seeing, 
which are complicated to model, potentially involving a large number of degrees of freedom. We 
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Fig. 6. — Distributions of ef , ejj , and cr^ Mb from 50 realizations. Left: 2100 stars are used for 
each realization. Right: 100 stars are used for each realization. The arrows point at the quadrature 
means of the distributions which are listed in TableEJ Note that the horizontal scales are different 
in the two panels. 



have seen in the SNAP case that adding 6 jitter degrees of freedom to the PSF model makes the 
PSF model errors twice as large as when we fit for only 5 misalignment parameters, so it seems 
likely that the PSF modeling performance will be degraded by larger number of parameters in the 
system. 
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